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Abstract 

An approach for understanding the behavior of multiplicity distributions in 
restricted phase-space intervals derived on the basis of global observables is pro- 
posed. We obtain a unifying connection between local multiparticle clusters and 
the scale-invariant power-law behavior of normalized factorial moments. The 
model can be used to describe multiparticle processes in terms of a decompo- 
sition of the observed intermittent signal into contributions from clusters with 
varying number of particles. 
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1 Introduction 



An inverse power-law dependence of normalized factorial moments (NFMs) F q (5) on 
the size of a phase-space bin 5 (intermittency phenomenon) obtained from multipar- 
ticle production in high-energy experiments can serve as a signal for dynamical local 
fluctuations with self-similar structure Q. For a single bin, this power-law behavior 
can be written as 

( n li}\ /An* 3 

F q (5) = oc ( T J , n W = n(n-l)...(n- 9 + l), (1) 

where n is the number of particles in the phase-space bin of size 5 in which these local 
multiplicity fluctuations are investigated. The angular brackets imply an averaging 
over all events in the sample and A is a full phase-space interval defined in rapidity, 
azimuthal angle, transverse momentum or a combination of these variables. The in- 
termittency exponents 4> q are related to the so-called anomalous fractal dimensions d q 
and the Renyi dimensions D q as follows 

d g = -^j, D q = D(l-d q ), (2) 

where D is the topological dimension of phase space. If particles are randomly dis- 
tributed in phase space, then (p q = and D q = D. For monofractal multiplicity 
distributions, d q = const. For a multifractal behavior, d q is a function of the rank q of 
the NFMs. 

The scale-invariant fluctuations manifest themselves as "spikes" - clustering of many 
particles in small phase-space bins for single events - and are common to many areas 
of physics (see for reviews). 

The factorial-moment method and the concept of intermittency in high-energy 
physics have been borrowed from the theory of turbulence. There, continuous prob- 
ability densities are used as the mathematical tool for a theoretical description of 
fluctuations |J. The main problem in such an approach is the comparison with exper- 
iments having a finite multiplicity of particles for single events. However, if statistical 
noise has a Poissonian nature, then the method of NFMs follows immediately due to 
its suppression of the statistical fluctuations caused by the finiteness of the number of 
particles in restricted phase-space bins. Then, the values of the NFMs are equal to the 
usual normalized moments obtained from probability densities. That is why, from a 
theoretical point of view, attempts have been made to understand intermittency also 
via the analysis of usual moments obtained from models borrowed from hydrodynam- 
ics |1]'IH31- However, little attention has so far been devoted to a systematic treatment 
of dynamical models involving a finite number of particles in phase space. 

In this paper, we develop the method for calculation of the NFMs making use of 
discrete probability distributions, i.e., a priori taking into account the finiteness of the 
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number of particles in a phase-space interval. Such an approach allows us to derive a 
multiplicity distribution in small phase-space bins from global characteristics of samples 
with correlations between particles. In this context, this paper may be considered as 
a continuation of the study of the projection method for the transition in probability 
distributions from full to restricted phase-space intervals HQ. 

The method is based on a natural assumption about the existence of multiparticle 
clusters in phase space and, for a given cluster-size distribution, can lead to an exact 
solution for the rise of NFMs with decreasing 5. Usually, some degree of arbitrariness 
exists in the definition of clusters. Here, a cluster in an individual event is considered in 
a general and traditional sense - as a bunch of many particles with a very small extension 
in the phase space under investigation. Let us note that such concept of spikes is wider 
than that of clusters, since occurrence of the cluster is caused by dynamical reasons 
while the spike may have a purely statistical nature. 



2 General Formalism 

First, let us define notations which will be used throughout this paper. The generating 
function G(z) for a multiplicity distribution P n and the unnormalized factorial moments 
F q are defined as follows 

oo 

G{z) = £ P n z n , F q EE (n®) EE G<*\z) | z =l • (3) 

n=0 

The normalization condition J2^=oPn = 1 leads to G(z = 1) = 1. Then, the NFMs F q 
for the multiplicity distribution P n are given by the relation 

* - 1 w 

A capital letter N will be used to specify the number of particles in the full phase space 
of size A, and n represents the number of particles in a restricted phase-space interval 
of size 5, so that n < N for 5 < A. 

Let us consider a collision between two particles yielding exactly iV final particles 
for each event in some full phase space of size A with a topological dimension D. Let 
us divide the full phase-space volume into M D non-overlapping bins of size 

s = w- (5) 

In this paper for simplicity we will consider the case of a flat phase-space distribu- 
tion, i.e., each infinitely small cell of phase space is regarded equally probable. In this 
case, none of the quantities characterizing fluctuations depend on the position of the 
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phase-space bin under study. To compare the model predictions with the experimental 
non-flat distributions, therefore, the transformation of a given non-flat distribution into 
a flat one should be performed [jnj . 

The intermittent fluctuations in physical systems may manifest themselves as lo- 
calized dynamical spikes (multiparticle clusters) in individual events. Let us mention 
briefly two well-known extreme cases of phase-space distribution: 

i) If all N particles are equally distributed in a given phase-space volume, then 



the multiplicity distribution for particles in 5 is a positive-binomial one with a 
generating function Gn(8, z) of the form 

G N (5,z) = (pz + g) N , g = l-p, (6) 

p = M- D = i (7) 
From (|]), we have the following form of the NFMs 

Nli] 



^ = — > g = l,2,...,iV. 



The extreme case thus gives <p q = d q = 0, D q = D. This means that all spikes observed 
in single events have a purely statistical nature. Note that F q < 1 for all q. However, 
for N —>■ oo, the distribution (|6]) tends to a Poissonian one with F q = 1. 

ii) Another (unlikely) situation occurs if all N particles group in one single point-like 
cluster for all possible events. Then the generating function is || 

G N (5,z)=pz N + g, g = l-p, (9) 

with p of the form (^). This distribution law emphasizes that, for a given bin size S, 
only two possibilities can occur: either all N particles are found in the bin or none. 
For the NFMs, one gets 

ATM /Ax 9 " 1 



f < = -mW ' « = ^--- N - ( 10 ) 

This is a maximum possible intermittency (4> q = d g = 1, D q = 0). From a geometric 
point of view, this case corresponds to a point-like object having topological dimension 
zero. 

The two examples presented above lead to the two main theoretical questions: how 
do such dynamical spikes lead to the actual form of intermittency ([!]) with a small (but 
non-zero) exponent 4> q and how does one construct the multiplicity distribution for this 
case. The aim of the present analysis is to derive a probabilistic scheme of the realistic 
intermediate situation with < <p q < q — 1 in ([!]), using known global characteristics 
of the sample. 
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In our approach we restrict ourselves to the following cases: 

1) All single particles (monomers) and multiparticle clusters are placed in phase 
space independently, without any correlations. Within the analytical model to be 
discussed below, all possible correlations are those between particles inside each cluster. 

2) We shall consider the simplified case of treating the clusters as "point-like ob- 
jects". In fact, we suppose that the probability that the cluster is emitted on the 
boundary of the restricted phase-space interval is very small, i.e., 5 ^> 5 min ~ O, 
where Q is the cluster size. Since intermittency is defined via the scale-invariant ratio 
A/5 = M D , we can rewrite this condition as 

£~-^->M* (11) 

Hence, for a given M, a cluster can be considered as a point-like object if 1) the size Q of 
the cluster is very small; 2) the total phase-space volume A of multiparticle production 
is large. For real experiments, in fact, condition ( |TTD means that we consider the case 
when the number of phase-space bins M is not very large. 

Consider the following structure of a fluctuation. Let us assume that not only single 
particles but also point-like multiparticle clusters can exist in single events. Let each 
of them contains exactly m particles (m-particle clusters), where m = 2, 3, 4 . . .. Then, 
only the two cases are possible: a cluster is inside 5 or outside of this interval. We 
assume that the cluster phase-space distribution is flat. If only multiparticle clusters 
with a fixed m exist, then the multiplicity distribution to have n particles (belonging 
to clusters) inside S has the binomial-like form: 



0, if n m s, s = 1, 2 . . . 

Cff p™g Nm ~m , otherwise, 



P n (m,N m ) = { ^ a „ . ' (12) 



where n — 0, 1, 2 . . . mN m , p is given by (0) and N m is the total number of m-particle 
clusters in full phase space. Cjy are the binomial coefficients. The N m is a constant 
for the sample of events and will be weighted over all possible samples below. The 



generating function for (|12"D is 



G { n}^ z) = £ z-C* N yg N — = (pz m + g)»~ (13) 

4 = 

As mentioned in the introduction, we will consider the case of no correlations be- 
tween different clusters and monomers. Then, the generating function for the prob- 
ability of having n particles in 8, if m-particle clusters (m = 2, 3 . . .) and monomers 
(m = 1) exist, is the product of the generating functions ( |i~3"D 

TV 

G R (5,z,N)= II (14) 

m=l 
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where Gjfi is the generating function (|6|) for uncorrelated monomers. Note that the 
generating function (|i4|) satisfies the normalization condition Gr(5, z = 1, N) = 1 by 
construction. The set R = {N m } = (Ni,N 2 , A^v) represents a cluster configuration 
(cluster-size distribution) which satisfies the following constraint 

N 

N=Y,™N m . (15) 

m=l 

Note that the extreme case Nn = 1, when all N particles form a single cluster- 
monomer, is possible only if N\ = N2 = . . . iV/v-i = 0. 

A real situation is expected to be more complicated, because, for a fixed N, the 
events can differ one from another by the cluster configuration R. In this case, R is a 
random variable. Then, one averages the Gn(5, z) over all possible configurations, i.e., 

00 

0(5, z,N) = J2 P c R on{ GB{5, z, N) = (G R (5, z, N)) R , (16) 

R=l 

where P^° nt is a probability distribution depending on iV and controlling the relative 
weight of the multiplicity distributions for different cluster configurations. 

Moreover, for actual multihadronic systems produced in high-energy experiments, 
one needs to average over all events with different N. Then, 

00 

G(6, z)=J2 PnG(5, z, N) = (G R (S, z, N)) RjN . (17) 

N=0 

Here, P/v is the known multiplicity distribution for full phase space. 

Thus, within the framework of the considered model, we have found the multi- 
plicity distribution in restricted bins from the following global observables containing 
dynamical information on the fluctuations: 

1) The multiplicity distribution P R on{ for the cluster configurations R. 

2) The multiplicity distribution P/v for particles in the full phase space A. 

Of course, both quantities have to be defined on the basis of physical knowledge on 
the dynamical structure of the particle clustering. 

According to (H) and (Hf), the NFMs for flTTD have the following form 



The quantities inside the angular brackets define (usual) factorial moments for a fixed 
cluster configuration R and a fixed total number of particles in the full phase space: 
q qi QN-2 

p (X\ - V V V C qi C q2 (7W-ip{l> ^{2} p{N-l} p{N} ,-,q\ 

r q\ )— 2-i ■ ■ ■ ^q ^qi ■ ■ - ^q N -2 r q-q\ r qi-q2 - ■ - r qN-2-qN-i r qN-i- 

<31=0<J2=0 <Jjv-i=0 
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Here, we have introduced the following definitions of the moments obtained from 
G£ } (M: 

JtW = iVf y, (20) 

F g {m} = EM^LpV^- m > 1. (21) 

It is convenient to introduce, instead of N m , the probability W m of any particle 
(chosen at random) to belong to an m-particle cluster. This probability can be written 

as 

JV 

W m = mN m /N, £W OT (iV) = l. (22) 

m=l 

The set of W m for m = 1, . . . , N forms the probability distribution which we shall 
call the cluster configuration distribution. Hence, R used above represents a particular 
form of the cluster configuration distribution W m . 

In fact, the quantities W m are identical to those used in the theory of percolation 
clusters (see JHJ for review). There, the cluster is defined as a group of occupied 
lattice sites connected by nearest-neighbor distances. Then N m /N can be treated as 
the average number (per lattice site) of m-particle clusters and W m defined from 
is the probability of any lattice site to belong to a m-particle cluster. 



3 Examples 

Though the theoretical attempt to describe a sample in terms of discrete distributions 
is still in a very early stage, one may already interested in its possible observable 
consequences. 

To illustrate the intermittent properties of the model, let us consider a specific 
system in which the particle multiplicity iV is fixed for all events and the P^ nf has a 
sharp maximum near R ~ R, i.e., 

^r="w={J;2?l (23) 

Then, 

G(6,z,N)z> Y[G { ^(6,z). (24) 

m=l 

where N m is the most probable cluster configuration. For a given N, expression 
permits us to calculate the NFMs if we know the form of N m (or the corresponding 
W m ). Then, from (E¥), we have 

FJ5) = (25) 
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where F q (S) is given by (|19D, if we substitute N m instead of N m . Note that ( |25| ) does not 
depend on the global particle density N/A, since, for any value of N, all dependence 
of the NFMs on the bin size is determined by the ratio A/5 = M D , but not by A itself. 

To be specific, we present the results of our calculation of the NFMs for a few 
particular cases: 

1 ) Contribution from two-particle clusters 

To study this case, we use the following simple cluster configuration 

W l = 0.4, W 2 = 0.6. (26) 
In Fig. HI we present the results of our calculation for F 2 to F 5 with N = 50 and 



D = 1 (closed symbols). For this value of N, according to (|22|), configuration (p6|) 
approximately corresponds to Ni = 20, N 2 = 15. 

From Fig. [I] we can see that the NFMs have a tendency to grow linearly with 
increasing of the number of bins M. The increase can be fitted by the power-law (fj) 
with a intermittency exponent <f) q and a coefficient of proportionality a q , 

In F q = <j) q \nM + \na q , (j) q > 0. (27) 

Since the main aim of this section is only to illustrate a real possibility of the inter- 
mittent rise of the NFMs with increasing M, we do not pursue the purpose to fit the 
results obtained by (^7]) for this (hypothetical) cluster configuration. 

2) Contribution from three-particle clusters 

The open symbols in Fig. p] represent the behavior of F 2 to F 5 for a sample with 
three-particle clusters (D=l). Here, 

Wi = 0.4, W 3 = 0.6. (28) 

This corresponds to Ni = 20, A^ 3 = 10 (for iV = 50). The sample with three-particle 
clusters has stronger intermittent behavior than that obtained using (p6|). 

4 Bunching-Parameter Method 



In this section, we argue and numerically show that the bunching parameters [13,14 



are more sensitive to the features of a cluster configuration than the NFMs discussed 



above. In fact, the result of this subsection complements the statement made in [I^j, 
where it has been shown that the NFMs are insensitive to the structure of fluctuations, 
i.e. rather different behaviors of the underlying local multiplicity distribution P n (5) 
with decreasing 5 can have the same intermittent trend of the NFMs. 
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The definition of the bunching parameters r] q is [p"3|.[T 



_ q P q (S)P q . 2 (S) 

where P n (5) is the probability of having n particles inside a restricted phase-space bin 
5. For example, for a Poisson distribution, one has rj q (5) = F q (5) = 1, so that qth order 
BP measures the deviation of the local multiplicity distribution from the Poisson. In 
general, independent particle production leads to T] q (8) ^ 1 for 5 — > 0. In this case, 
however, the bunching parameters are ^-independent constants. For example, if all 
spikes are purely statistical for events with a fixed finite multiplicity N, i.e., when the 
cluster configuration reduces to the trivial case 

W m = 6 mA , (30) 

then ([14]) becomes a positive binomial distribution and the corresponding BPs have 
the following ^-independent form 

q-l-N 

*< = T=2=N' (31) 

The interest in the use of BPs to extract information on the cluster configuration 
lies in the fact that the BP of rank q is sensitive only to multiparticle clusters with 
m < q particles. That is, a given BP acts as a filter for clusters having a small number 
of particles. This follows directly from the definition (p9|). On the contrary, the NFM of 
rank q is sensitive to spikes with n > q particles and acts as a filter for clusters with large 
number of particles. This means that the bunching parameters have a complementary 
property which is very important to obtain a refined insight into various multiparticle 
systems with intermittent behavior. 

To demonstrate this point, we shall use (|4]) and the relation P q (5) = G <yq \z) \ z=0 
I q\. Note that for this calculation we can use the expression (|I9~D, substituting P g (m, N m ) 
(see (|l2f) ) multiplied by the factor q\, instead of the Fj m >. 

In Fig. |2| (a),(b), we present the behavior of the NFMs for W\ = 0.4, W 3 = 0.6 
(only monomers and three-particle clusters exist) and W\ = 0.22, W 2 = 0.4, W 3 = 0.3, 
W4 = 0.08. The latter configuration represents a sample with two-, tree- and four- 
particle clusters. The behavior of the NFMs is the same for these rather different 
samples, not only qualitatively, but also quantitatively (actually, these configurations 
were chosen for an illustrative purpose, because they exhibit the same behavior of the 
NFMs). In contrast to the NFMs, the behavior of the BPs is different for the same 
cluster configurations (see Fig. |3] (a) (b)). This means that the bunching parameters 
provide a simple and effective method to analyze models and compare them with 
experimental data. 



This point has already been illustrated in |14[ using a Monte-Carlo model, where we 
have shown that systems with a similar trend of the NFMs exhibit different behavior 
when one studies them with the help of bunching parameters. 
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5 Conclusion 



We have presented a statistical analysis of multiparticle fluctuations in ever smaller 
phase-space bins for different dimensions using the characteristics of primarily observ- 
able - multiparticle clusters. This new theoretical method links the global characteris- 
tics (a cluster configuration distribution and a multiplicity distribution of particles in 
the full phase space) of the sample to local multiplicity fluctuations of a multiparticle 
process. As we have seen, the theoretical approach developed in its simplest form (no 
correlations between clusters) reveals many promising features. The model can lead to 
an intermittent rise of the NFMs with decreasing 5, even for the simplest special cases 
of the global characteristics taken for illustrative purposes. In our model, we "a pri- 
ori" take into account the clusters, without any physical examination of the clustering 
phenomenon itself, which, up to now, is a standard topic in soft hadron physics where 
the application of QCD is very restricted. 

One of the reasons for the approach presented in this paper is to give an exact 
method to derive the intermittent behavior in various cluster models. In particular, 
the general method developed here can be used to study intermittent behavior of the 
statistical models |0| of a fragmentation process in usual phase-space variables, rather 



then for mass distributions. 

There is an interesting point connected with our approach. In fact, our prescription 
is to decompose the local fluctuations into clusters with varying number of particles. 
Apparently this situation is closely analogous to that in which the intermittent signal 
is decomposed into contributions from clusters living on specific scales (the so-called 
method of orthogonal wavelet transform used in high-energy physics ||16|| ). Note, how- 
ever, that the connection between that method and our approach is not trivial, since 
the conceptual distinction between these strategies is mathematically clear-cut. 
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1 2 3 4 5 

In M 

Figure 1: NFMs as a function of the number of bins M . Solid lines represent the 
configuration W\ = 0.4, W 2 = 0.6, dashed lines correspond to W\ = 0.4, W 3 = 0.6. For 
both cases N = 50. 
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Figure 3: Dependence of In rj q (M) (q = 2, . . . , 5) on In M for the same cluster configu- 
rations as in Fig. [| 
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